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Abstract. The lazy random walk, where the walker has some probability of staying 
put, is a useful tool in classical algorithms. We propose a quantum analogue, the 
lackadaisical quantum walk, where each vertex is given I self-loops, and we investigate 
its effects on Grover’s algorithm when formulated as search for a marked vertex on the 
complete graph of N vertices. For the discrete-time quantum walk using the phase flip 
coin, adding a self-loop to each vertex boosts the success probability from 1/2 to 1. 
Additional self-loops, however, decrease the success probability. Using instead the 
Shenvi, Kempe, and Whaley (2003) coin, adding self-loops simply slows down the 
search. These coins also differ in that the first is faster than classical when I scales less 
than N, while the second requires that I scale less than Finally, continuous-time 
quantum walks differ from both of these discrete-time examples—the self-loops make 
no difference at all. These behaviors generalize to multiple marked vertices. 


PACS numbers: 03.67.Ac 

1. Introduction 

Random walks, or Markov chains, are the basis for a variety of classical algorithms [T]. 
One typically wants the random walker to move readily, but there are some scenarios 
when it is desirable for the walker to be lazy, meaning it has some probability of staying 
put. For example, if a normal (non-lazy) random walker starts in one of the two vertex 
sets of a bipartite graph, then at each time step, the walker will only be in one vertex 
set. By making the walk lazy, however, the walker can have probability of being in 
both vertex sets, improving its coverage of the graph. Such lazy random walks equate 
to adding self-loops to the vertices of the graph with appropriate weight, and they have 
been utilized in a variety of classical algorithms, including PageRank |2], graph covering 
|3], and image processing |1]. 

Given the success of lazy random walks in the classical regime, we propose in this 
paper a quantum analogue called lackadaisical quantum walks, dehned to be a standard 
quantum walk [SI El El on a graph with the addition of I self-loops to each vertex of the 
graph. So the greater I is, the more the walker prefers to stay put. Note this differs 
from the “lazy” quantum walk proposed by |8], hence our choice of a different name. It 
also generalizes three-state lazy quantum walks on the line 0 EHl El E2], which only 
have one self-loop at each vertex. Since there exist both discrete- and continuous-time 
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Figure 1. (a) The complete graph with N = 6 vertices. A vertex is marked, as 
indicated by a double circle. Identically evolving vertices are identically colored and 
labeled, and the labels indicate the subspace basis vectors that the vertices belong to. 
(b) With I self-loops at each vertex. 


varieties of quantum walks , we will consider both in our investigation of lackadaisical 
quantum walks. 

In particular, we explore the addition of self-loops to one of the best-known 
problems in computing: unstructured search, whose famous quantum solution is 
Grover’s algorithm To review, given a “database” with entries 1, 2,..., iV, and 

an oracle f{x) that outputs 1 for a particular “marked” entry w and 0 otherwise, a 
classical computer expects to query the oracle N/2 = Q(N) times before finding w, 
since it could be the first guess or the last. In the quantum setting, |1), |2),..., |A^) 
are computational basis states, and the oracle R^, = (—acts by flipping the phase 
of a marked basis state |tc) while leaving the rest unchanged, i.e., Rwlw) = — |w) and 
Rw\x) = |a;), Vx 7^ w. This reflection through |w) can be written as R^ = I — 2\w){w\. 
Grover’s algorithm ra finds |ta) in only Q{'/N) queries, which is a quadratic speedup 
over the classical Q{N), and it does so by initializing the system in an equal superposition 

1 ^ 

is) = ^Ei^) 

over the basis states and repeatedly applying 

I^ c-L < 


( 1 ) 


where = 2|s)(s| — I is a reflection through Applying these two reflections 

tt'/N /4 times, the state is rotated from |s) to |w) with probability near 1. 

As a graph problem, we formulate this unstructured search problem as search on 
the complete graph of N vertices for a particular marked vertex, an example of which is 
shown in figure la Since each vertex is connected to every other, there is no structure 
demanding an order to which we visit the vertices. Thus a classical random walk that 
jumps from vertex to vertex, checking at each step if it has found w, expects to make 
Q{N) such steps and checks to hnd the marked vertex. A quantum walk, on the other 
hand, searches in Grover’s Q{'/N), and we explicitly show this in the next two sections 
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in both discrete and continuous time na. Then we add I self-loops to each vertex, as 
shown in hgure lb, and show how the lackadaisical walk affects the search algorithms. 
For the discrete-time quantum walk, we get an improvement in the success probability 
of the algorithm with a particular “coin” with / = 1. Additional self-loops, or search 
with a different coin, hurts the algorithm. The continuous-time algorithm, on the other 
hand, is not affected by the self-loops at all. Finally, we generalize the results to multiple 
marked vertices. 


2. Grover’s Algorithm as a Discrete-Time Quantum Walk 

We begin by analyzing quantum walk search on the complete graph with no self-loops, 
in discrete-time in this section and continuous-time in the next. For both discrete- and 
continuous-time quantum walks, the quantum walker jumps from vertex to vertex, and 
the vertices of the graph label computational basis states of an A-dimensional “vertex” 
Hilbert space C^. For discrete-time quantum walks, however, this is insufficient to dehne 
a local unitary operator [151 [IB], so we necessarily include an additional d-dimensional 
“coin” Hilbert space supported by the directions/edges that the particle can jump 
along from each vertex. For the complete graph, each vertex is connected to the other 
N — 1 vertices, so d = N — 1. Let |s^) and |sc) be equal superpositions over each space: 

N d 

1=1 ^ 1=1 

Then the system |'^) begins in the equal superposition over the entire ® Hilbert 
space: 

IV’o) = ki-) ® kc)- 

The quantum walk is dehned by repeated applications of 
Uo = S-{In® Co) , 

where Cq is the “Grover diffusion” coin im 
Co = 2kc)('Sc| ~ Id^ 

and S is the “flip-flop” shift [T8| that causes a particle to hop and then turn around, 
e.g., a particle at vertex 1 that points towards vertex 2 will, after an application of S', 
be at vertex 2 and point towards vertex 1: S(|l) 0 |1 —)■ 2)) = |2) 0 |2 —)• 1). Note that 
l-^o) is the equilibrium distribution of this walk, so Uo\'ipo) = I'ko)- 

To turn this into a search algorithm, we apply a different coin Ci to the marked 
vertex and still use Co on the unmarked vertices ra. so the search operator is 

U = S ■ [(/at - |tc)(tc|) 0 Go |w)(tc| 0 Cl]. (2) 

Two choices for Gi are common [18]. The hrst is Gf^'’ = — Go, which causes U to become 

f/ = S- [(Jjv-2|u;)(u;|)0Go] 
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= S ■ {In® Co) ■ {In - 2\w){w\) ®Id 

'-V-' '-V-' 

Uq Rm 

= Uo-{R^®Id). (3) 


Note Rw is the phase flip in Grover’s algorithm, so this applies a phase flip to the marked 
vertex followed by a step of the quantum walk. The second common choice for Ci is 
^SKW _ which was hrst introduced by Shenvi, Kempe, and Whaley dl to solve 

search on the hypercube, and later used by Ambainis, Kempe, and Rivosh [TH] to solve 
search on arbitrary dimensional periodic square lattices. 

With either of these coins, the system evolves such that there are only two types of 
vertices, as shown in hgure [T^ In particular, the marked red a vertex evolves differently 
from the identically-evolving unmarked white b vertices. Since the a vertex can only 
point towards b vertices, and the b vertices can either point towards the a vertex or 
other b vertices, the system evolves in a 3D subspace, and we take equal superpositions 
of these vertices/directions as the basis vectors: 

\ab) = \a) ® ^ |a ^ 

^ S ®\b^a) 


In this {\ab), \ba), \bb)} basis, the initial state is 

1 


iV’o) = 


Vn 


ab) + \ba) + \/N — 2\bb) 


In general, the search operator (j^ is different for and But for the complete 

graph, they are identical; with either coin, the search operator is 

/ 0 — cos 9 sin 9 \ 

(4) 


U = 


-1 

0 


0 0 
sin 9 cos 9 


where 


cos 9 = 


N 

N -V 


and sin 9 = 


2^/]V^ 
N-I 


Repeatedly applying this operator to the initial state, the success probability evolves as 
shown in hgure for N = 1024 and 2048 vertices. We see that the success probability 
reaches 1/2, and the runtime (he., number of applications of U to reach the maximum 
success probability) scales less than linear (he., classical) since doubling N results in a 
runtime that is less than double. In particular, we expect it to scale as 0(a//V) to be a 
formulation of Grover’s algorithm. 
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Figure 2. Success probability as a function of the number of discrete-time 

applications of U (with either the or coins) for search on the complete 

graph with N = 1024 (solid black) and 2048 (dashed red) vertices. 


To prove this behavior and hnd the precise runtime, we hrst hnd the eigenvectors 
and eigenvalues of U. They are 


IV’+) = 


1^-) = 


1 -b cos 9 
3 + cos 6 


1 -b cos 9 
3 -b cos 9 


1 

2V1+COS 9 


COS 

9 

— 

iV3 

+ 

cos 

9 


1 

2V1+COS 9 

[y/l — 

cos 

9 

+ 

iy/3 

+ 

cos 




1 







/ 

1 

2\/l+cos 9 

1_a/1 — 

cos 

9 

+ 

iVS 

+ 

cos 

9 

\ 

1 

2V1+COS 9 

1_a/1 — 

cos 

9 

— 

iy/3 

+ 

cos 



E~ = 


, 


E_ = e-'^^ 


IV’-i) = 


1 — cos 9 
3 -b cos 9 


/_ Ji+co^ \ 

1 —COS d 


1 +CQS^ 

1 —COS 9 


, E_i — —1 




where (j) is defined such that 

, 1 -b cos 9 


cos 0 = 


, and sin 0 = 


a/( 1 — cos 9) (3 -b cos 9) 


2 ' " 2 
Now we express the initial state in terms of the eigenvectors of U. Consider 


1 , , , ,, 1 /I -b cos^ 

^ “ 71V3 + COS0 


/ /1—cos 6 \ 

' l+cos6 


1—cos 6 


1+cos 6 

\ 2 

For large iV, sin^ 2/y/N implies that 9 ^ 2/y/N, so the first two components of this 
are 

1 -COS0 [¥J2 19^ 9 1 


4 2 


1 -b cos 9 


2 
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which means the last term dominates for large N. That is, 

Since l-^o) ~ \bb), the system after t applications of U is 

U‘m ~ T {u‘ \i>+) + U‘ |v>-)) = ^ (e'*‘ \i>+) + e-f’ |«i_» . 

When (pt = 7 i/ 2, i.e., 



TT 


TT 


TT 


2 sin 


-1 


yQ ( 1 —cos6»)(3+cos6») \ 2 sin ^(0 /a/2) 2a/2 




the state of the system is approximately 


1 

71 


(* IV’+) 


1^-)) 


1 /1 + cos 6 

a/ 2 V 3 + cos 6* 


/ 

V 


3+cqs 6 
1+COS0 

/ 3+cos^ 
' 1+COS0 


0 


\ 


I 

71 


/1\ 

-1 

\V 


So the system roughly evolves from \bb) to being half in \ab) and half in \ba), which 
from the \ab) component gives a success probability of 1/2. This agrees with figure]^ 
the success probability reaches 1/2 after 7rA/l024/2A/2 36 and 7rA/2048/2A/2 50 

applications of U. Repeating the algorithm an expected constant number of times to 
boost the success probability near 1, the algorithm still finds the marked vertex in 
0(a/]V) applications of t/, which is the same scaling as Grover’s algorithm. 

This explicit proof that the success probability reaches 1/2 in 7r\//V/2-\/2 
applications of U seems original, even though it is well-known that the success 
probability reaches 1 in continuous-time [H] (which we review next), or in discrete¬ 
time with one self-loop per vertex [18] (which we review later). Note that the first 
discrete-time quantum walk search algorithm Cl was search on the hypercube, not the 
complete graph. 


3. Grover’s Algorithm as a Continuous-Time Quantum Walk 


In continuous-time, quantum walk search on the complete graph with no self-loops 
(Z = 0) was previously investigated by [Tl|, and we review the results here. Continuous¬ 
time quantum walks do not require the coin space, so they walk in the iV-dimensional 
Hilbert space supported by the vertices of the graph. The system begins in the 
equal superposition over the vertices: 


1^(0)) = Is,;), 

and evolves by Schrodinger’s equation 

with Hamiltonian 


H = —jA — \w){w 
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Figure 3. Success probability as a function of time for continuous-time search on 
the complete graph with N = 1024 (solid black) and 2048 (dashed red) vertices, at the 
critical 7 = 1/N. 


where 7 is the jumping rate (he., amplitude per time), A is the adjacency matrix of the 
graph {Aij = 1 if i and j are adjacent and 0 otherwise), and |tn) is the marked vertex 
we are looking for. The hrst term effects a quantum walk [TT] while the second term 
acts as an oracle [19], with 7 setting their relative strength. 


As shown in hgure la, there are only two types of vertices: the marked vertex and 
the unmarked vertices. So we take equal superpositions of them to be basis vectors of 
a 2D subspace: 


|a) 

\h) 


red) 


1 


E i»)- 

zGwhite 


In this {|a), |&)} basis, the initial state is 


IV'(O)) = |s^) = i|a) + \ 




N 


and the Hamiltonian is 

H = -'j 


N-2 


(5) 


When 7 takes its critical value of 1/N [HI [201 1^ - evolving by Schrodinger’s equation 
with this Hamiltonian yields the success probability shown in hgure [^ We see that it 
reaches a maximum value of 1 , and the runtime (which should be Grover’s ©(vF)) 
scales better than linear (classical) since doubling N less than doubles the runtime. 

To show this analytically, note that the eigenvectors of the Hamiltonian with 
7 = 1 /Af are \i/j±) cx |si,) =F |a) with corresponding eigenvalues E± = — 1 ± 1/a/F +1/N. 
Then the initial state is |'0(O)) = (iV’-i-) + \i ’-))and the marked vertex is |a) = 
(— |'0+) + |'0-)) /a/2- Since the Hamiltonian is time-independent, solving Schrodinger’s 
equation yields 


|^(f))=e-^^|^(0)) 
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= e 


-\Ht 


72 


^ ( 7 +) + \' 4 ’-)) 


= |^+) + 7 -)) , 

where AE = — E_. When AEt = tt, i.e., 

TT TT 


t 

AE 2 ’ 


this becomes 


\^jJ(t)) = e (- |?/^+) + IV'.)) = e 


_ —i_E_7r/AE 


a). 


Thus the system evolves to the marked vertex with probability 1 in time t:\/N /2, which 
agrees with hgure with N = 1024 and 2048; as expected, the success probability 
reaches 1 at time 7r7l024/2 k, 50.265 and vr72048/2 71.086. 


4. Discrete-Time Quantum Walk With Self-Loops 


Now we include I > 0 self-loops at each vertex, as shown in figure lb As before, there 
are only two types of vertices, a and b. But now the a vertex can also point towards 
itself, so the system evolves in a 4D subspace spanned by equal superpositions of the 
vertices / directions: 


\ab) = \a) ® 

^ ^ ® I& ^ a) 

In this {|aa), \ab), \ba), \bb)} basis, the initial state is 

I'/'o) = — , (\/l\aa) -|- VN — l\ab) + VN — l|6a) 

7Ar(iV + /-l) V 

+ 7(Ar-l)(Ar + /-2)|66)). 


With self-loops, the = —Co and CSKw ^ coins now result in different search 
operators (|^ and evolutions, which we analyze separately. 

With = —Co, [15 showed that with this coin and I = 1 self-loop at each 
vertex, two applications of the search operator ([^ corresponds exactly to Grover’s 
iterate Q on the vertex space. Reproducing their argument, the vertex and coin 
spaces have equal dimension N in this case. Then Co = Rgj- = 2|s)(s| — In is 



















Grover Search with Lackadaisical Quantum Walks 


9 



Figure 4. Success probability as a function of the number of discrete-time 
applications of U with the coin for search on the complete graph with N vertices 
and I self-loops at each vertex. The solid black, dashed red, and dotted green curves 
correspond to iV = 1024 with Z = 1, 2, and 3, respectively, and the dot-dashed blue 
curve corresponds to = 2048 with 1 = 2. 


the reflection about the equal superposition |s) over the iV-dimensional computational 
basis in Grover’s algorithm. With this substitution, the search operator (|^ becomes 
U = S ■ {In ® -Rs-l) ■ {Rw ® In)- Acting by this on the initial equal superposition state 
iV’o) = k) ® |s), we get = S{Ryj\s) ® Rs±\s)) = ® Rw\s)- Acting a second 

time, U'^l'ipo) = S{Rii;Rg±\s) ® Rs±Rw\s)) = Rs±Rw\s) ® RwRgQs), which is precisely 
Grover’s iterate Q on the first tensor factor. Thus the success probability reaches 1 in 
TTy/N /2 applications of f/, which is an improvement over the success probability of 1/2 
without any self-loops. 

To And the behavior of the algorithm for general I > 0, consider the search operator 
Q or ([^, which in the {|aa), |a6), \ba), \bb)} basis is 


U = 


( cos^ 
0 


• sin d 
0 


— sin Q — cos Q 


\ 


0 


0 


0 

— cos ( 
0 

sin d) 


where Q is defined such that 

N-l-\ 

cos u = 


\ 


0 
sin 
0 

cos (j) 


N + l-V 
and d> is defined such that 


and sin^ = 


2V/(iV-l) 

N + l-1 ' 


( 6 ) 


cos 0 = 


N + l-3 
N + l-V 


and sin 0 


2^N + 1-2 
N + l-1 


Repeatedly applying this to the initial state, the success probability evolves as shown 
in figure 1^ Beyond / = 1, additional self-loops causes the buildup of success probability 
to stall, resulting in a lower maximum success probability. The figure also reveals that 
the maximum success probability only depends on I, which is reasonable since the I = 0 
and I = 1 cases achieve success probabilities of 1/2 and 1, independent of N. 
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To find the precise behavior of the algorithm, we hnd the eigenvalues and 
eigenvectors of the search operator U. The eigenvalues are 


-l,l,e--,e'“, 


where 


cos a = 


cos 9 + cos (j) N — 2 
2 ^ N + l-V 

+ cos 6 * + cos 0) (2 — cos 6 * — cos 0) a/ {2N + / — 3)(/ + 1) 


iV + / - 1 


and the corresponding (unnormalized) eigenvectors are 


sin^ l + cos^ 1 + COS0 l + cos0 

1 + cos 9 sin 0 ’ sin 0 ’ sin 0 ’ J 


1 + cos 9 sin 0 sin 0 sin 0 
sin 9 1 + cos 0 ’ 1 + cos 0 ’ 1 + cos ( 


sin 6 * cos 6 * —cos 0 .sina cos 6 ^ —cos 0 .sina ^ 
sin 0 ’ 2 sin 0 sin 0 ’ 2 sin 0 sin 0 ’ 


/ sin 6 * cos 9 — cos 0 . sin a cos 9 — cos 0 . sin a 

= I ^ T) 7 ) ^ 7 * ^ T) 7 ^ 7 1" * ^ 7) ^ 

V sm 0 2 sm 0 sm 0 2 sm 0 sm 0 


Then 


^ 1 + cos 0 sin^ 6 * J ^ J 

2:j—-7 . 2 < 

1 + COS sm 0 

..^1 —cosa ..^1 —COSO; 1 —cosq;\^ 

= (0)2 ^ - ,2 ^ - ,4- — j 

0 sm 0 sm 0 1 — cos 0 / 

is dominated by the last term because sin 0 Ri 0 in the denominator of the second and 
third terms is small for large N, which implies that 1 — cos0 ~ 0^/2 in denominator of 
the last term. Thus if we normalize it to leading-order, 

1 — COS0 / 1cos 0 sin^0 \ ^ ^ i,,\ 

77^-X 2—-+ 0+a + 0-a ~ (0, 0, 0, 1)^ = |66). 

4(1 —cos a) \ 1 + cost^sm 0 / 

Note that the initial state |0o) ~ 1^^)- Then after t applications of U, the system is in 
the state 


(0,0,0,1)T = |66). 


l^o) 


1 — cos 0 / 1 + cos 0 sin^ 9 


4(1 —cos a) \ l-|-cos 6 *si 


>1 e^'^* 0 -|_Q. -\- e 


We choose t such that at = n, i.e., the runtime is 

t = — ^ < .rr/om-l f VW+‘^) \ J _ 

a 


n/ sin ^ 


I = cN 
I = uj{N) 
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for large N. At this runtime, the state of the system is 

1 — cos d) 


l^o) 


4(1 — cos a) 
1 — cos (j) 
4(1 — cos a) 


''H-a 


„ 1 + cos (j) sin^ 9 , 

2 -— - K - 01 

1 + cos 6 sin^ 0 

sin 6* 2 — 3 cos 6 + cos 0 

_ 4- -- 


sm( 


sm( 


2 — 3 cos 6 + cos ~ cos 6 + cos < 

I ; 1 ^" 


sm (p 1 — cos ( 

Then the success probability p is given by the sum of the squares of the hrst two terms: 
(1 — cos 0)^ /16 sin^ 6 * + (2 — 3 cos 9 + cos 0)^' 


P = 


16(1 


cos a 


sm 


Plugging in for cos a, sin0, cos0, sin 6 *, and cos 6 *, this is 

I = o{N) 


P = 


16/(A^- 1) + (3/- 1)^ 
4(/ + l)2(Ar + /-2) 


(/+ 1)2 
16+9c 1 
4c(c+l) N 
9_ 

4/ 


/ = cN 


I = u:{N) 

for large N. These expressions for t and p agree with £gure|^ for search with N = 1024 
vertices and I = 1,2,3, the runtimes are respectively 7rA/l024/ \/2{l + 1) 50,41,36 

with corresponding success probabilities 4//(/ + l)^ = 1, 0.889, 0.75. With N = 2048 and 
I = 2, the runtime is 7rA/2048/■\/2(2 + 1) 58 with success probability 4(2)/(2 + 1)^ 


These results indicate that the maximum success probability decreases as I 
increases. Even so, there is still an improvement over the loopless success probability of 
1/2, so long as / < 5. For I beyond this, the success probability is less than 1/2, but if I 
is any constant, the overall runtime is still Grover’s When I scales greater than 

a constant, the quadratic quantum speedup is lost due to the classical repetitions of 
the algorithm needed to boost the success probability. Despite this, when I = o{N), we 
still obtain an improvement over the classical algorithm’s Q{N) runtime. Finally, when 
I = Ll{N), the success probability fails to increase beyond its initial scaling of 0(1/A^), 
and so the quantum algorithm is no better than the classical one. 

Now consider = —Id- The search operator ([^ in the {|aa), \ab), \ba), \bb)} 

basis is 


/-I 

0 

0 

0 ^ 

0 

0 

— cos 9 

sin^ 

0 

-1 

0 

0 

yo 

0 

sin 9 

cos 6 * y 


(7) 


where 


cos 9 


and sin 9 

N + l-V 


2VN + 1-2 

N + l-1 


Clearly, |aa) is an eigenvector of U with eigenvalue —1. The remaining part of U 
corresponding to |a 6 ), | 6 a), and \bb) takes the same form as U for the I = 0 case (|^. 
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Figure 5. Success probability as a function of the number of discrete-time 
applications of U with the coin for search on the complete graph with N = 1024 

vertices and I self-loops at each vertex. The solid black, dashed red, dotted green, and 
dot-dashed blue curves correspond to Z = 0, y/N = 32, 2N = 2048, and = 32768, 
respectively. 


Since 6 is small for large N, those results carry over: the success probability reaches 1/2 
in 

TT 

t = - y 

1—cos &) (3+cos &) 

l sm I ^ 


applications of U. The scaling of this with N depends on the scaling of Z. In particular, 
for large N, 

{ l = o(N) 

[ 55 , 7 / 

An example of this evolution is shown in hgure with the success probability reaching 
the expected 1/2 at time 7r\/1024/2A/2 36 for both I = 0 and I = y/N = 32, time 

7r^rT2Vl024/2y2 62 for I = 2N = 2048, and time 7rV32768/2y2 201 for 

I = Ar3/2 ^ 327g3_ 

With this coin, self-loops cause the success probability to take longer to reach 1/2. 
As long as the number of self loops scales less than or equal to N {i.e., I = 0{N)), the 
runtime still scales as Grover’s 0(a//V). Furthermore, there is still a speedup over the 
classical algorithm so long as I scales less than (he., I = o{N^)). Compared to the 
Cf’*’ coin, which obtains Grover’s Q{y/N) runtime when I = 0(1) and a speedup over 
classical so long as I = o{N), this coin is in some sense more robust to self-loops; 

it takes more of them for the walk to lose its quantum speedup. 


5. Continuous-Time Quantum Walk With Self-Loops 

Let us see how the continuous-time quantum walk algorithm is affected by the presence 
of I self-loops at each vertex. If we count each self-loop to contribute 1 to the diagonal 
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of the full iV-dimensional adjacency matrix so that An = I at vertex i [H then in the 
two-dimensional subspace spanned by {|a), |&)}, the Hamiltonian is 


hr = -7 


N + l-2 j 


Note this is simply the Hamiltonian with no self-loops ([^ plus II. Adding a multiple 
of the identity matrix to the Hamiltonian in this manner constitutes a rezeroing of 
energy or an overall phase, so it has no observable effects. Thus the self-loops do not 
change the evolution at all; with or without self-loops, at the critical 7 = 1/N, the 
success probability reaches 1 at time iry/N /2. Thus the continuous-time quantum walk 
algorithm is completely robust to lackadaisical errors in our model using I self-loops at 
each vertex. 


6. Generalization to Multiple Marked Vertices 


All of these results are straightforward to generalize to the case of k marked vertices. 
We assume that k = o{N) since the number of marked vertices cannot scale more than 
the number of vertices, and ii k = cN, then one can classically find a marked vertex 
in a constant number of guesses. The classical search would take an expected Q{N/k) 
time to find one of the k marked vertices on the complete graph of N vertices. As for 
the quantum walk, let us consider each of the cases above. 

Beginning with discrete-time quantum walks, with k > 1 marked vertices 
and I self-loops at each vertex, the system evolves in a 4D subspace spanned by 
{|aa), \ab), \ba), \bb)}. With the Cf''’ = —Co coin, the search operator ([^ or (|^ is 

0 \ 
sin 0 
0 

cos (f) y 


^ cos 9 — sin 0 0 

^ 0 0 — cos 0 

— sin 6 — cos 6 0 

0 0 sin 0 

where 9 is defined such that 


cos 9 


N -2k-I+ l 
N + 1-1 


and 


sin 9 


2^{N -k){k + l-l) 
N + 1-1 


and 0 is defined such that 


cos 0 = 


N -2k + l-l 
N + l-1 '■ 


and 


sinffl = 


2^yk{N - k + l - 1 ) 
N + l-1 ' 


This has the same form as the case of one marked vertex (§ , and since 0 is small for 
large N, the solutions carry over: define a such that 

cos 0-1-cos 0 N — 2k 


cos a = 


N + l 


f Note that some treatments count each self-loop as 2 in the adjacency matrix, which would result in 
An = 21, but it makes no difference to our result. 














Grover Search with Lackadaisical Quantum Walks 


14 



Figure 6. For search on the complete graph with N = 1024 vertices, A: = 16 marked 
vertices, and I self-loops at each vertex, the success probability as a function of the 
number of discrete-time applications of U with the coin. The solid black and 
dashed red curves are respectively I = 1 and 32. 


Sin a = 


a/( 2 -b cos 9 + cos 0) (2 — cos 9 — cos 0) 

2 

y/{2N -2k+ 1- l){2k + 1-1) 

A^ + / - 1 ■ 


Then after 


TT 

t = — 
a 


-^/N 


, . ( a/c(c+2) 

TT/sm 


I = o{N) 

I = cN 
l=uj{N) 


applications of U, the success probability reaches a maximum value of 
(1 — cos 0)^ /16 sin^ 6* -b (2 — 3 cos 9 + cos 0)^' 

p = 


16(1 —cos \ sin^ 0 

k [16N{k + l-l) + 9{l- If - 4k{l - 1) - I2k‘^] 


4{2k + l-iy{N-k + l-l) 

I = o{N) 
l = cN 
I = uj{N) 

1 ^ where the success probability reaches 
13 and 4 ■ 16(16 + 32 - l)/(2 ■ 16 + 32 - l)^ s 
~ 9 applications of U, as expected. 


(2k+l-ip 
16-I-9C 
4c(c-|-l) N 
9k 
41 

figure 


This is shown 
7rVl024/^2(2-16 + 1-1 ^ 
at 7r^I(M/^2(2-16 + 32-l 

With the = —1^ coin, the search operator ([^ in this basis is 


1 at 
0.758 


/-I 

0 

0 

0 ^ 

0 

0 

— cos 9 

sin^ 

0 

-1 

0 

0 


0 

sin 9 

cos6' y 
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Figure 7. Success probability as a function of the number of discrete-time 
applications of U with the coin for search on the complete graph with N = 1024 
vertices, fc = 16 marked vertices, and I self-loops at each vertex. The solid black and 
dashed red curves are respectively I = 1 and 32. 


where 


cos 6 


N -2k+ 1-1 
N + l-1 


and sin 6 


2^k{N - k + l-1) 
N + l-1 


This has the same form as the case of one marked vertex Q, and since 9 is small for 
large N, the solntions carry over: we reach a snccess probability of 1/2 in 


t = 


TT 


2 sin 


-1 


—cos 6)(3-l-cos 9) 
2 


2+2k 

TTy/c+l 


y/iV 


2+2k 


2+2k 


VI 


o{N) 

cN 

uj{N) 


applications of U. This is shown in figure where the success probability reaches 1/2 at 
7r-\/l024/2-\/2 ■ 16 ~ 9 and 71^/2 + l-\/l024/2-\/2 ■ 16 ~ 15 applications of f/, as expected. 
Finally for the continuous-time quantum walk, the Hamiltonian is 

H = -J + \ 

\ -s/kiN -k) N-k + l-1 I' 


This simply adds II to the Hamiltonian with no self-loops 1221 , which is a rezeroing of 
energy or global phase, so it has no observable effects. 


7. Conclusion 

We have proposed a quantum analogue of classical lazy random walks, called 
lackadaisical quantum walks, where the quantum walker has some preference to stay 
put by introducing I self-loops at each vertex of the graph. We have investigated the 
consequences of this for quantum search on the complete graph, showing the self-loops 
can have vastly different effects on quantum search depending on the type of quantum 
walk. For discrete-time quantum walks with the Cf''’ coin, one self-loop per vertex 
boosts the success probability, but additional self-loops hinders it. With the coin. 
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however, rather than hindering the maximum success probability, the self-loops slow 
down its buildup. This hindrance is less potent, eliminating all speedup over classical 
search when I = u{N‘^) compared to the hrst coin’s I = u{N). Continuous-time quantum 
walks, on the other hand, are not affected at all by the self-loops. All of these results 
extend to multiple marked vertices. 
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